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Abstract 

The inclusion process is a stochastic lattice gas, which is a natu- 
ral bosonic counterpart of the well-studied exclusion process and has 
strong connections to models of heat conduction and applications in 
population genetics. Like the zero-range process, due to attractive 
interaction between the particles, the inclusion process can exhibit a 
condensation transition. In this paper we present first rigorous results 
on the dynamics of the condensate formation for this class of models. 
We study the symmetric inclusion process on a finite set S with total 
number of particles N in the regime of strong interaction, i.e. with 
independent diffusion rate m = mjsi — > 0. For the case Ntun — > oo we 
show that on the time scale 1/mjv condensates emerge from general 
homogeneous initial conditions, and we precisely characterize their lim- 
iting dynamics. In the simplest case of two sites or a fully connected 
underlying random walk kernel, there is a single condensate hopping 
over S as a continuous-time random walk. In the non fully connected 
case several condensates can coexist and exchange mass via interme- 
diate sites in an interesting coarsening process, which consists of a 
mixture of a diffusive motion and a jump process, until a single con- 
densate is formed. Our result is based on a general two-scale form 
of the generator, with a fast-scale neutral Wright-Fisher diffusion and 
a slow-scale deterministic motion. The motion of the condensates is 
described in terms of the generator of the deterministic motion and 
the harmonic projection corresponding to the absorbing states of the 
Wright-Fisher diffusion. 
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1 Introduction 



The inclusion process, introduced in pQ as the dual of a model of heat 
conduction, and further developed in [3J, is a natural counterpart of the 
well-known and extensively studied exclusion process. The inclusion pro- 
cess on a finite set S combines random walk jumps at rate mp(i,j)/2 with 
inclusion jumps, where each particle at site i is attracted by each particle 
at j independently with rate p(i,j). Here p(i,j),i,j £ S are the rates of 
a symmetric irreducible continuous-time random walk on S. The inclusion 
jumps introduce a form of attractive interaction between the particles. 

When m tends to zero and simultaneously many particles are in the 
system, due to the attractive interaction between the particles, large piles 
of particles will be formed at individual sites. In fact we proved in [3] 
condensation for the stationary measure, in a limit where both m — >• and 
the number of particles N — >• oo. In this paper we prove dynamical results 
for the condensation phenomenon. More precisely, we study how clusters 
arise from an arbitrary initial condition, and how they move and merge into 
a single condensate, which then jumps over the finite lattice as a random 
walk. We therefore make m = rriM dependent on N € N such that itin — > 
as N — > oo, which implies that diffusion will get slower and the attractive 
inclusion interaction will create condensates. 

This model can be interpreted as a multi-allele version of the Moran 
model [3], describing the evolutionary competition of several species in a 
fixed size population. The inclusion part describes reproduction and death 
and m jv plays the role of an additional mutation rate, which is typically very 
small on the reproduction time scale (see e.g. [5] and references therein). 
It is also important as a dual process to models of energy or momentum 
transport [fl [2j El [7] . 

Related recent theoretical results include explosive condensation in a to- 
tally asymmetric model [8j which exhibits a slinky motion of the condensate 
also observed recently in [9J for non-Markovian zero-range dynamics. Re- 
sults on the equilibration dynamics in zero-range processes with decreasing 
rates are currently under investigation [TU]. In order to have a well de- 
fined limit dynamics we require symmetry of the p(i,j), but our proof is 
not based on potential theoretic methods as used in [111 [T2] for zero-range 
processes. Due to the system-size dependence of the diffusion rate we can 
apply more general multi-scale methods (see e.g. p3]). In particular, we 
use a two-scale decomposition of the generator and construct the generator 
of the limit process similarly to results in p3] . 

The decomposition consists of a Wright-Fisher type diffusion part which 
runs at "infinite speed" (in the limit N —¥ oo), and a p(i, j')-dependent drift 
part. The limiting motion is then described by the harmonic projection of 
the drift part on the absorbing set of the diffusion. In the simplest case 
of two sites, or similarly, in the fully connected case where all p(i,j) are 
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strictly positive, the absorbing set of the diffusion is the set of corner points 
of the simplex E = {x € [0, l] s : Ylies Xi = -*■}• This corresponds to the 
immediate formation of a single condensate, which then hops over the set S 
as a random walk. In that case we can characterize the limit dynamics in 
a relatively simply way as the solution to a martingale problem for linear 
functions. In the more general case, in a first stage, several condensates will 
form, and interact via intermediate sites. Condensates do not split but can 
merge, and this coarsening dynamics eventually leads to a single condensate 
which then again moves over S as a random walk. In this case, we compute 
the generator explicitly making use of the separation of time scales and the 
martingales for the fast Wright-Fisher diffusion. 

The rest of our paper is organized as follows. In Section [2] we introduce 
the model and present the main results and possible extensions. In Section 
Owe collect general results on tightness and the comparison of semigroups 
and work out the simple two-site and fully connected cases. The proof of 
the general case is treated in Sections H] and Sections [6] where we give a 
general definition of the limiting generator of a two-scale system in terms of 
projection operators on harmonic functions. In Section[5]we give an explicit 
computation of the limiting generator in the case of 3 lattice sites, which 
already contains most of the difficulties of the general case. 

2 Definitions and main results 

2.1 The model and the auxiliary slow-fast system 

We consider a finite set S, with associated jump rates p(i,j) = p(j, i), j,i 6 S 
of a symmetric, irreducible random walk. 

The symmetric inclusion process (w(t) : t > 0J with parameter m > 
based on p(i,j) is then defined as in [2] to be the continuous-time Markov 
process on the configuration space Z+, with generator 

*$ IP M = £ PiiJH (~ + m) {f(v ij ) - m) , (i) 

ijes 

where r\i G Z + denotes the number of particles at site i € S, and rf£ = 
Vk ~ $i,k + k S S denotes the configuration where one particle moved 
from i to j. Since the total number of particles is conserved, we consider the 
process with generator ([T]) on the state space Qn = {v £ ^+ : J2ies = 
with a fixed total number N G N of particles. 

We consider m = mjv dependent on the number of particles N such that 
mjv — > as N — >• oo, which implies that diffusion becomes slower and the 
attractive inclusion interaction will create condensates. We are interested 
in the limiting dynamics of the rescaled process after accelerating time by a 
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factor 0jy ;= a/rajy with a > 0. We will assume throughout the paper that 

niN — > , NrriN — > oo so that — ^- = — — >• . (2) 

N m^N 

We denote by ei,i G 5 the canonical unit vectors of h(S) with entries 
e i(j) = <5y> Consider the rescaled process (x N (t) : t > 0) defined by 

a^(i) := {rii(0 N t)/N:i€S). 

This is a Markov process on the simplex E = jx G [0, l] s : X]jg5 x « = 1} 
with generator 

^v/(x) = 8NP{hj)Nx t (f + iVx,) ^f( x -j r e i +j r e j ) - /(x)) . (3) 
yes 

Assuming smooth / in ([3]), Taylor expansion of the right-hand side gives, 
using the symmetry oip(i,j), 

■&nS(x) = - ^2 jP(i,j)(xi ~ Xj)(dijf)(x) 
yes 

+ \ E p{h3)^e N {dy){x) + 0(6 N /N) , (4) 

i,jeS 

where we abbreviated 

Here, by ([2]), the correction term denoted by 0(6n/N) is a function of x 
which converges to zero uniformly in x, as N — > oo. We will show that 
this correction terms can be ignored and one can study the easier auxiliary 
process (y N (t) : t > 0) on E with generator 

L]\f := L + 0nL' 

= ~ ~ + pi^^yiVjdfj . (6) 

yes yeS 
As iV -> oo this is a slow-fast system with two scales, where the fast part 

£' = ^Ep(m)kv 3 4 (7) 

yes 

corresponds to a Wright-Fisher diffusion with absorbing states given by 

.c/ := {x G £ : p(i,j)xiXj = for all i,j G S 1 } . (8) 
The slow part is a deterministic motion with generator 

L = ~Y1 f J0(W - • ( 9 ) 

yes 
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The limiting dynamics of the process with generator Ln is intuitively de- 
scribed as follows. The fast (#ArL')-part of the generator pushes "infinitely" 
fast to the set of absorbing states whereas the slow (L) part tries to move 
away from s/. So the motion consists of infinitesimally moving away from 
$rf and immediately being projected back onto stf '. This will result in an 
"effective, projected process" on the set s£ ' . We will describe this procedure 
of "projection" rigorously and in more detail in section FOl below. 

2.2 Main results 

Our main results rigorously characterize the limiting dynamics of the process 
(x N (t) : t > 0), by showing convergence to a limit process (x(t) : t > 0) 
which concentrates on the subset srf C E. Furthermore, the set of corner 
points 

<T := {a : i € S} C $4 (10) 

is in turn absorbing for the limit process, as we can conclude from the first 
theorem. In the following we speak of weak convergence on path space, if for 
every T > the processes (x N (t), < t < T) converge to (x(t), < t < T) 
weakly in the Skorokhod topology. 

Theorem 2.1. 1. Assume that x N '(0) — »■ x° € in distribution. Then 
(x it) ■ t > 0) converges weakly on path space to (x(t) : t > 0), which 
is a Markov chain on ^ with initial condition x(0) = x° and generator 

Af(e l ) = ^^Mf^)-f(^) ■ (11) 

jes 

2. In the fully connected case, i.e., if p(i,j) > for all i,j E S, then the 
same holds for a general initial conditions. I.e., if x N (0) — > x° € E 
in distribution, then (x (t) : t > 0) converges weakly on path space to 
(x(t) : t > 0), the process on with generator pi[) . and with initial 
condition P(x(0) = a) = E[x°]. 

So if the process is started from a configuration in i.e., when (as 
N — > oo) all the mass concentrates on a single site, all future configurations 
are of this type and the single pile of mass (condensate) performs a random 
walk on S with rates proportional to p(i,j). The same holds in the fully 
connected case for general initial conditions, since in that case srf = 
and only the corner points are absorbing for the Wright-Fisher part. For 
general symmetric p(i,j), has a more complicated structure and the limit 
process is not ergodic on si '. In the following we will describe the dynamics 
for general initial conditions under the additional assumption 

p(i,j)£ {0,1} ioraRiJeS, (12) 
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i.e., where edges are either connected or not with uniform jump rate 1. 
The more general case can be treated as well but is more complicated to 
formulate, and we will comment on that in a discussion in Section 12.31 Let 

p(i, j) = (1 - p(i,j)) ^ P(i, k)p(k, j) > (13) 
fees 

be the number of two-step connections between sites i, j £ S, which are not 
directly connected. 

To state the result we also need to introduce the harmonic measure v x 
on s/ for the Wright-Fisher diffusion with generator L' given in (JJJ), with 
x £ E as initial condition. This is the distribution of the point x(oo) € £/ 
where the diffusion with generator V is absorbed when started from x £ E 
(see Section r4.2l below for more details). For the initial condition in Theorem 
12.11 item 2, note that for non-random x° the harmonic measure is exactly 
u x°( e i) = x i> since the marginal dynamics on each site is given by a diffusion 
(see Section 13.31 for details) and only the corner points are absorbing. For 
the general case below there is no easy explicit formula for u x o. 

Theorem 2.2. Assume IWi) and assume that at time zero x N (0) — > x° £ 
E in distribution. Then (x N (t) : t > 0) converges weakly on path space 
to (x(t) : t > 0), which is a Markov process on with initial condition 
x(0) ~ E^o] an d generator 

jcs \ies ) \ ies / 

The generator of this process consists of two parts. The first part 
corresponds to an effective Wright-Fisher diffusion between sites which are 
connected by a path of length two only, which results from exchanging mass 
via the intermediate site. This mechanism of mass exchange is also studied 
heuristically in [8] for a totally asymmetric model. The second part corre- 
sponds to all the neighbouring mass of an empty site accumulating on that 
site with a rate proportional to the total neighbouring mass. In the simplest 
case where there is only one neighbour with non-zero mass, this corresponds 
to a jump of a single pile. Therefore piles can exchange mass continuously 
or merge, but they never split and the set of corner points ^ where all mass 
concentrates in a single pile is absorbing. For the inclusion process this ab- 
sorbing set is reached on the same time scale 6n as the stationary motion on 
the corner set takes place. This is different from other models such as the 
zero-range process, where these dynamics happen on different time scales 

nans mi 1121129. 
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Figure 1: Simulation illustrating the main results by a sample path of the 
process with normalized occupation numbers. The lattice is S = {1, ... ,4} 
with nearest neighbour connections and periodic boundary conditions. Ini- 
tially, there is diffusive mass exchange between sites which are not directly 
connected, and later all mass concentrates on a single pile. Fluctuations are 
due to finite parameter values, where N = 1000, to at = 0.01 and On = 100. 

Note that the initial condition for the limit process is not x° and con- 
vergence in the usual sense does not hold at time 0, since the process is 
"immediately" projected onto the absorbing set si '. The initial condition is 
then distributed as the absorption point on &/ starting from x°, which is 
consistent with the right limit x(0+). 

2.3 Discussion 

The general limit dynamics of Theorem 12.21 is illustrated in Figure [2] for a 
ring of 4 sites. After initial diffusive mass exchange between sites 1 and 
3, the system reaches the absorbing set of corner points and turns into a 
jump process. Figure [2] illustrates the fully connected case (second item 
in Theorem 12. ip where the system instantaneously reaches a corner and a 
single pile performs a random walk on the lattice. 

Our results hold for symmetric jump rates p{i,j), since otherwise the 
Taylor expansion @ would contain first order drift terms diverging with N, 
which would lead to infinite limiting speed of the clusters. Therefore the 
limit dynamics are not well defined unless the p(i,j) are symmetric. 

In the formulation of Theorem l2.2l we restrict ourselves to p(i, j) 6 {0, 1}, 
only for simplicity of presentation. Our methods can be applied straight- 
forwardly to the general case, only the concrete computations in Section [5] 



7 




Figure 2: Simulation illustrating the main results for a fully connected lattice 
S = {1,2,3} of three sites, where immediately all mass concentrates on a 
single pile. The left panel shows the evolution of all sites and the right 
panel the evolution of site 1 for longer times. Fluctuations are smaller than 
in Fig. ffl since here N = 10000, m N = 0.001 and 6 N = 1000. 



have to be adapted. Assuming p(l,2) = p(2, 1) = p, p(2, 3) = p(3, 2) = q 
and p(l,3) = p(3, 1) = for the case of 3 sites, the continuous part of the 
limit generator A (|14p then takes the form 

^(q-p)x\X3 di3 + ~(qxi+(p-q)xl) (px 3 +(q-p)xl) df 3 . 

Compared to (|66|) and (|14|) , there is an additional drift term and the diffusive 
term is no longer of Wright-Fisher type. This can be formulated for general 
lattices analogously. The jump part of A still takes the same form as in (|14p 
on a general lattice. 

3 A technical lemma and first results for the two- 
site and the fully connected case 

This section consists of two parts. 

1. In section 3.1, we show the following. If the auxiliary process (y N (t) : 
t > 0) with generator Ljy ([6]) converges to a limiting process (y(t) : 
t > 0) and if the "true process" (x N (t) : t > 0) with generator «2V 
([3]) converges along a subsequence, then this limit is the same process 
(y(t) : t > 0). As a consequence, if the sequence (x (t) : t > 0) is 
tight, then it has the limit characterized by the auxiliary process. This 
observation will be proved in Lemma 13. II below and allows us to work 
with the simpler generator Ln to characterize the limit process. 

2. In subsection 13.21 below we first consider the case of two sites, which 
is easier than the general case because the limiting process is a pure 
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jump process on the set of unit vectors {e±,e2}- Then we work out 
the fully connected case, which is similar to the two site case, and 
where the limiting process is also a pure jump process on the set of 
unit vectors {e^ : i € S}. 



3.1 General results on comparison and tightness 



We start with a technical lemma that we will apply later to justify dropping 
the O(0n/N) part in the Taylor expansion of the generator 

Lemma 3.1. Suppose J>?n and Ln are generators of Markov processes on the 
compact state space E, such that on a common core of functions o€ C C(E) 
we have 

||^v/-^/||oo<|||/|||/37v. (15) 

Here ||.||oo denotes the supremum norm, /3tv is a sequence of numbers con- 
verging to zero as N — > oo, and ||.|| denotes a norm such that for all f S o€ 
and t > 



sup sup | 

N 0<s<t 



< OO . 



(16) 



Then we have: 

1. For the associated semigroups, for all i > and f 6 

lim \\e tLN f - e^/Hoo = . 

N^foo 



(17) 



2. IfYmx N ^ 00 e tLN f = S t f for all t > 0, f G C(E) with (S t : t > 0) a 



Feller semigroup, then lim^r- 



Nf = s t f as well. 



3. If the processes with generators J£n and Ln are tight on the path space 
-D(0,oo) of rcll functions (right- continuous with left limits), then they 
have the same limit process. 

Proof. Pick / £ r Jtf . Start from integration by parts (or Dyson's) formula 
(see e.g. [17] p. 367) and use that e sJ ^ N is a contraction in the supremum 
norm to estimate 



f 



/lie 



< " '(JSfjv - L N )e {t - s)LN f ds 



< 



< 



|| (Sf N - L N )e^ LN f Woods 



P. 



,(t-s)L N 



f\\\ds 



< fi^t sup I r 

0<s<t 



(18) 



as N — > oo for all t > 0. 
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For the second item, let / 6 C(E), and for given e > choose g G 
such that ||g — /||co < £• Then for N large enough we have, by item 
1, \\e tIjN g — e N g\\ 00 < e and by the contraction property of e tLN ,e N , 
l|e* L ^ (/-<?) Hoc, ||e^(/-5)||oo <e. This implies \\e tL " f - e^/IU < 3e 
and if one of the terms converges to a semigroup, so does the other. 

Finally, the third item follows because a Markov process on the path 
space -D(0, oo) is uniquely determined by its semigroup. □ 



Lemma 3.2. The sequence of Markov processes (x (t) : t > 0) on a com- 
plete separable metric space E is tight on the path space Z?(0, oo) of rcll 
functions if 

1. The sequence [x (t) : t > 0) is stochastically bounded on D[0,oo), 
and 

2. For all f £ C£°(E), denoting compact smooth functions (which are 
uniformly dense in C(E)) and all e > there exist non-negative ran- 
dom variables Z]\f(5,f,e) such that 

E[f(x N (t + u)) - f(x N (t))\x N (t) =x)\< E[Z N (6,f,e)\x N (t) = x] 

(19) 

with probability 1 for all t > e and < u < 5, and 

lim lim sup E[Z N (S,f,e)] = . (20) 

The conditions are sufficient but not necessary, for a proof, see e.g. [TU] . 
Lemma 3.11 or [20], Chapter 3. The first condition will always be ful- 
filled since our state space E is a compact, finite dimensional simplex and 
1 1^(2) I |oo < 1- The second condition has been modified for our purposes, 
and ensures boundedness of the modulus of continuity of sample paths start- 
ing the process from any strictly positive time and will be verified later. At 
time the sequence of processes we consider is actually not tight, since in 
the limit it gets instantly projected onto an absorbing state of the Wright- 
Fisher diffusion L' . So we can show convergence only for all positive times 
t > 0, and by right-continuity of sample paths we can extend the limit pro- 
cess uniquely to all t > 0. We will show in Section 0] that this extension is 
consistent with the harmonic measure, describing the absorption of the fast 
Wright-Fisher diffusion. 

3.2 The two site case 

We start now with the simplest case 

S = {1,2} and p(l,2)=p(2,l) = l, (21) 
where S contains only two sites and £? = = {e±, e2}- 
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Theorem 3.1. Assume l21\) with initial conditions x N (0) such thatx N (0) — » 
x° € [0, 1] in distribution as N — > oo. Then the process (x N (t) : t > 0) con- 
verges weakly on path space to (x (t) : t > 0) , which is a simple random walk 
on {ei,e2} with jump rates a/2 and initial condition x(0) G {ei,e2} with 
P(x(0) = ei) = x\. 

The rest of this section is devoted to the proof. Since S consists of two 
elements we abbreviate x\ = x, and because in the process x\ + X2 = 1 is 
conserved, we effectively have only one variable (denoted by x). The Taylor 
expansion leading to then gives in terms of x: 

^Nf(x) = (-±9 N m N (2x - \)d x + 9 N x{l - x)d'i) f{x) + O(0 N /N) 

= : L N f(x) + 0(9 N /N) (22) 

where the terms 0{9^/N) go to zero uniformly in x, as N — > oo. The 
generator 

L N = -\a{2x-l)d x + 26 N x{l-x)d 2 x (23) 

appearing in (|22p is a Wright-Fisher diffusion with mutation where in front 
of the diffusion term x(l — x)d x we have a factor 9n, i.e., the diffusion part 
(also called "genetic drift part" in the population dynamics language) of the 
process is accelerated w.r.t. the mutation process. 

This implies intuitively that the process will tend to evolve "immedi- 
ately" towards the fixed points of the genetic drift process, which are the 
homo-zygotes x = 1 and x = 0, and the mutation process will then lead to 
a random flipping between these two states. 

First, we show that we can apply Lemma 13. II to our case. Since this is 
slightly technical, we state it separate lemma. 

Lemma 3.3. The generators and Ln of (|22[) satisfy the assumptions 
of Lemma \3. 1\ with \\\f\\\ = (||/"||oo + ||/'"||oo) with the set of poly- 
nomials as a common core. 

Proof. The approximate generator reads (for simplicity we put a/2 = 1 
in the following) 

L N = -(2s - l)d x + 2 6 N x(l - x)d 2 x 
and for the difference with the true generator we have the estimate 

\\s?Nf - L N f \\oo < c^nriu + Eunice) < ^(iiriioo + iiriu) 

from a standard expression of the remainder in Taylor expansions. The 
norm |[|/||| = (||/"||oo + ||/ w ||oo) can thus be used in (fT5|) . 

So in order to see that (|16p holds for this choice of norm, it is suffi- 
cient to prove that for the auxiliary process y N (t) with generator L^, the 
expectations 

®yf(y N (s)) = e sLN f(y) 
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have bounded second and third derivatives in x, uniformly in < s < t and 
N € N, for sufficiently many /. Choosing for / the polynomial f n (y) = y n , 
we have, 

LnIh = 29 N n(n - l)(/ n -i - f n ) + n(f n -i - fn) ~ nfn 

This can be interpreted as follows: the first two terms form the genera- 
tor of a Markov chain (nt : t > 0) that jumps from n to n — 1 at rate 
2 8^n(n — 1) + n corresponding to the classical ancestral dual of Wright- 
Fisher diffusion (Kingman's coalescent + mutation contribution), and the 
last term is an extra killing part. For more details on that connection see e.g. 
the nice account on duality in population models [18J. Using the Feynman- 
Kac formula we then obtain, 

iKN,t,n,y) :=E y f n (y N (t))=E%{e-f>° ds f nt (y)) 

where denotes expectation for the process on N with generator 

KNg{n) = N n(n - l){g{n - 1) - g{n)) + n(g(n - 1) - g(n)) . 

As a consequence, since the process with generator Kn only lowers the 
starting n in the course of time, we obtain the uniform estimate 

d k 



—iP(N,t,n,y) = E^(e-/o« S V 1 -V(n t -l)...(n ( -Hl))|<n J 



for all y £ [0,1], which gives ()16p . Finally, it is clear that on the compact 
simplex E the polynomials form a common core of the generators Ltv,J2?at. 
□ 



So Lemma 13. II guarantees that we can identify the limit from the auxil- 
iary process (y N (t) : t > 0) with generator Ln- We first show that the limit 
process concentrates on the set {0, 1} and then characterize it as a solution 
to the martingale problem of the claimed limiting process. 

Lemma 3.4. We have for allt > 

sup E y [y N {t){l - y N {t))] ^0 as N -> oo , (24) 
ve[o,i] 

and furthermore 

limsup#7V sup E y [y N (t)(l -y N (t))} <C (25) 

7V->oo y£ [0,1] 

for some C > 0. The same holds for x N (t)(l — x N (t)) of the original process 
with generator ££n . 
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Proof. We abbreviate f(y) = 2y(l — y), g(y) = 2(2y — l) 2 and compute 

jE y f(y N (t)) = E y (g(y N (t))) -26 N E y f{y N (t)) . 

Putting ifj t (y) = E y f(y N (t)), (p t (y) = E x g(y N (t)), we rewrite this equation 
as 

— ip t = tp t -2 6 N ip t , 
at 

which has the standard form of a slow-fast system (see e.g. [13] for more 
background), omitting the fixed argument y. This can be rewritten as 

-28 N t pB N t 

i>t = 4>oe- 26Nt + J <p s /e N e 2s ds . 

Using now the fact that y £ [0, 1], we obtain that <p < 2 and ip < 1 which 
gives 

< Mv) < e~ dNt + j- (l - e- 2dNt ) (26) 

which goes to zero as N — > oo for all t > 0, thus proving (|24p . The estimate 
(|26|) also implies the second statement ([25]) . 

The same argument holds for x N (t)(l — x (£)) if we replace yj( by ^ + 
0{9n /N) including the remainder terms, since we only use boundedness of 
ipt for N large enough. □ 

This shows that the limit process (if it exists) concentrates on the state 
space {0, 1} for all t > 0, and the estimate can also be used to show tightness 
using the criterion of Lemma l3.2l which implies existence of a subsequential 
limit by Prohorov's theorem (see e.g. [20], Chapter 3). 



Lemma 3.5. The processes (y N (t) : t > 0) and (x N (t) : t > 0) are tight on 
the path space D(0,oo). 

Proof. As mentioned before, the state space E = [0, 1] is compact and the 
first condition of Lemma [3 . 2 1 holds clearly. For the second condition, we first 
note that the process is not tight for t = since it is immediately projected 
onto the absorbing set si which is {0, 1} in this simple case. For all t > 
we can show tightness for the auxiliary process first by noting that for any 
smooth function / on [0, 1] 

E(f(y N (t+u))-f(y N (t))\y N (t)=y) <e( f t+ \L N f{y N {s))\ds\y N {t)=y) , 

(27) 

and we use 

rt+S 

Z N (5J,e) = S np \L N f(y N (s))\ds (28) 

t>e Jt 
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in the notation of Lemma 13.21 Then 

L N f(y) = -(2y - l)/'(y) + 9 N y(l - y)f"(y) (29) 
so that for all 5 > 

/ r t+s \ 

Ey( J \L N f(y N (s))\ds) < SWfWoo + 5C||/"||oo , (30) 

where we can use the estimate (|25p since t > 0. For every fixed e > this 
bound holds uniformly in N, and uniformly in the initial condition y £ E. 
It tends to as 5 \ since / is smooth and all derivatives are bounded on 
[0,1]. This implies tightness by Lemma 13.21 Since e > was arbitrary we 
have tightness for all strictly positive times t > 0. 

The same argument holds for the original process (x N (t) : t > 0), since also 
third derivatives of / are bounded and additional terms on the right-hand 
side of d30]) are of the form 5 1 1 /"' || «, (C + N /N) . □ 

In order to identify the limit process by a martingale problem, the fol- 
lowing lemma gives a simplified version for a process that jumps with rate 
A back and forth between and 1. 

Lemma 3.6. Let Xt be a Markov process with values in {0, 1} and A > 0. 
If 

M t :=X t -X - I \(2X t - l)ds (31) 
Jo 

is a martingale, then Xt is the process that jumps with rate A between 1 and 
0. 

Proof. Denote by 

Lf(x) = Ax(/(0)-/(l)) + A(l - aO(/(l)-/(0)) = -A(2x - l)(/(l)-/(0)) 

the generator of the process jumping back and forth between 0, 1 at rate A. 
Let / : {0, 1} -+ M, then f(x) = xf(l) + (1 - x)f(0) = x(f(l) - /(0)) + /(0) 
is always linear. Therefore, using (|3ip 

f(X t ) - f(X ) - f Lf(X s ) = M t f (32) 
J o 

is a martingale for all / : {0, 1} — > R, which characterizes the process. □ 

Now we can conclude the proof of Theorem 13.11 Applying the generator 
Ln (|22p of the auxiliary process to the function f(y) = y gives that 

M N (t) := y N (t) - y N (0) - | jf* (2^( S ) - l) ds 
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is a martingale. Therefore, if (y N (t) : t > 0) — > (x(t) : t > 0) along a 
subsequence, for the limiting process we have that 

M(t) := x(t) - x(0) f (2x{s) - l) ds 

is a martingale. Thus the limit is the claimed jump process on {0, 1}. By 
Lemma 13.11 this is also the limit of the original process x N (t) . 

Since the limit process (x(t) : t > 0) has right-continuous paths the 
consistent initial condition is given by 

x(0) : = lim lim y N (e) , (33) 

which concentrates on {0, 1} by Lemma 13.41 To see that this has the dis- 
tribution claimed in Theorem 13.11 denote tp(t,N) = E y N(n\(y (t)). Then 
using the generator ([23]) . we see that 

*M1 = -(2,( ( , W ,- 1 ) 

which gives 

(p (t,N)=y N (0) + l(l-e- 2t ) . 

Taking the limits N — > oo and then t \, 0, and using that the limiting process 
is concentrated on {0, 1} gives the initial condition as claimed. 

Remark 3.1. We will see in Section^ that in general the limiting semi- 
group S(t) has the property that 5(0) = P is a projection onto harmonic 
functions on E. So if fi is the distribution of limiting initial condition x° as 
given in Theorems \ 2.1\ and \2.ty, the initial distribution of y(0) for the limit 
process is given by fiP = E^fz/^o] since (Pf)(x) = E^f/] for all f € C(E,M). 
For a symmetric diffusion on [0, 1] the harmonic measure v x on {0, 1} is then 
given by v x = x\, i.e. equal to the first component of x° £ E (see the second 
item of Theorem \2. 

3.3 The fully connected case 

We will see in the next section that Lemma f3.ll the concentration Lemma f3.4l 
and the tightness argument in Lemma [3.5l can be extended straightforwardly 
to the general case. In case the set S is fully connected, i.e. all p(i,j) > 0, 
we have s& = c € and the limit process still concentrates on the corner points 
and can be characterized by a simple martingale problem. This will provide 
a proof of the second statement of Theorem 12.11 

The following lemma gives the martingale characterization of the process 
jumping with rates p(i,j) from ej to ej. 
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Lemma 3.7. Let (y(t) : t > 0) be a process on {ei,i € S}, and denote 
m(t) = y(t).ei = t(y(t) = (%). If for all i 

M\t) := yi (t) - w (0) -Vf ( yj (s) - Vi(8))p(i,j)d8 (34) 

ieS Jo 

is a martingale, then y is the random walk jumping at rate p(i,j) = p(j,i) 
between ej and ej . 

Proof. As this is a straightforward extension of Lemma 13.61 we will be 
brief. On the corner points c <o of the finite simplex E a general function 

can again be written as a linear combination of simple linear functions yi, 
which implies the simpler martingale characterization along the same argu- 
ments as Lemma 13.61 in the two-site case. □ 

The generator Ljy = J2i jes © applied to the linear function fk(y) = 
Vk gives 

LNfk(y) = -^p(k, i)(xi - x k ) 
Therefore, for all k € S 

!/?(*)- v?(o) - £ f^M-y^) 

ies Jo 

is a martingale. So the limit along any convergent subsequence of (y N (t) : 
t > 0) will fulfill the condition in Lemma [3. 71 and is therefore unique. Again, 
the initial conditions of the right-continuous limit process can be defined as 
in (f33|) and have distribution Efz^o] by the same arguments, which proves 
the second statement of Theorem 12.11 

4 The general case 

4.1 Concentration and tightness 

Recall the generators ([J]) and © of the original and the auxiliary process, 

Ea 
-p[i,j){Xi - Xj)(dijf)(x) 

i,j€S 

+\ P{h3)^N{dlf){x) + 0(8 N /N) 
i,jes 

= Lf(x) + 6 N L'f(x) + O(0 N /N) . (36) 
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Individual terms for fixed i,j are precisely of the form dealt with in the two- 
site case. Therefore, Lemmas 13,41 and 13.51 on concentration and tightness are 
relatively straightforward to generalize, and we only give short proofs in the 
following. 

Lemma 4.1. We have for all t > for the auxiliary process 

p(i,j) supE^yf (t)yf (t)] ^ as N -> oo , (37) 

ydE 

and furthermore 

p(i,j) limsupfliv sup E y [yf(t)yf(t)} < C (38) 

AT->oo ye[0,l] 

for some C > 0. The same holds for the original process (x N (t) : t > 0) 
with generator J2?at. 

Proof. Fix i,j € S with p(i,j) = p(j, i) > (otherwise there is nothing to 
prove) . Analogous to the proof of Lemma 13.41 we can compute 

L'iviVj) = -2p(i,j)yiyj (39) 

and 

HviVj) = -j{yi'}2p(j, k )(yj -Vk) +Vj^2p{i,k)(yi -yiS) ■ ( 40 ) 
kes kes 

The second term is bounded in absolute value by the constant 



C := a max y~^(p(i, k) + p(j, k)) 



J k£S 

and only the first term is relevant, since it is multiplied by 6n- Using the 
notation 

<p t (y)=E y [y?(t) y y(t)] and Mv) = ®v [Hyf yf )] (41) 
we get the estimate 

< My) < e~ 2p(i ' j)eNt + , C — ( 1 - e -MiJ)6 N t\ (42) 

with the same arguments that lead to (|26p . The rest of the proof follows 
analogously to the two-site case. □ 



Lemma 4.2. The sequences [x N (€) : t > 0) for the full process and (y (t) : 
t > 0) for the auxiliary process are tight on the path space D(0, oo) of rcll 
functions. 
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Proof. We use again the criterion in Lemma 13.21 on the compact state 
space E and an analogous argument as in the two-site case. The only new 
element is the estimate of the generator analogous to (|30p . which follows 
from 

supE v [L N f(y N (t))} < V tejOII^/IU + CH^/lloo) (43) 

since all yf (t) G [0, 1] and for t > we can use the estimate ([38]) . Again, / is 
smooth on the compact space E so that all partial derivatives are bounded. 
□ 

So the limiting process exists for all positive times t > and concentrates 
on the absorbing set stf of the fast Wright-Fisher diffusion. The initial 
condition on s$ will be given consistently by right limits of the sample 
paths which coincide with the absorption probabilities of the Wright-Fisher 
diffusion, as we will see in the next subsections. 

4.2 Slow-fast Markovian systems and projection of Markov 
processes 

In the general case, the generator ([6]) still consists of a fast Wright-Fisher 
diffusion part and a slow mutation drift part. The limiting process will 
take place on the absorbing set s$ of the Wright-Fisher diffusion, which 
in the case of general p(i,j) is a richer set than the set of corner points 
= {ej : i G S}. 

We first define a the limit motion of a general two-scale Markovian sys- 
tem in greater generality, following closely the results in [TJ] . Let us consider 
(X(t) : t > 0), [X'(t) : t > 0) Markov processes on a compact metric space 
E with generators L, L' . Both generators should have a common core 
(which is given by the smooth functions C°°(E) in our case) and the corre- 
sponding semigroups e and e tL are defined on the Banach space C(E, R) 
of continuous functions with the supremum norm ||.||oo- Suppose X[ has 
a Borel measurable set of absorbing states g/ C E such that the hitting 
probabilities from any point x € E are well defined and given by 

u x (A) := lim F x (X'(t) € A) (44) 

t— >oo 

for measurable subsets Ads/. Of course we have v x = 5 X for all x G 

In the spirit of [T3], Theorem 1.11, we define the X' -harmonic projection 
operator 

P : C(E,R) -> Jt?(E) where Pf(x):=f f{a)v x (da) . (45) 

for all / G C(E, R), where the range of P is the set of L'-harmonic functions 
J^{E). Here, v x (da) is the harmonic measure for the generator L' on the 
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set si defined in (|44p . Hence, the function Pf is harmonic for L' and solves 
the Dirichlet problem Vip = with ip{a) = /(a) for all a G In particular 
Pf(a) = f(a) for a 6 <sz/, and as a consequence P 2 / = Pf, i.e., P acts as a 
projection from the set C(E,M) onto the set of harmonic functions. 

Denote by S K (t) the semigroup of the two-scale process X K with gener- 
ator L + kL'. The following result describes the limiting process as k — > oo 
by a semigroup on the harmonic functions. 

Theorem 4.1. Suppose that e tL ' f — > Pf as t — »■ oo /or a// / e C(P,]R) 
and i/ia£ £/ie operator 

A : C(E,R) -> iw/iere Af:=PLf (46) 

generates a Markov semigroup S(t) on Jff(E) (the range of P). Then 
Sn(t)f S(t)f as K — > oo /or a// / € Jf?(E) and t > 0. 77ie limiting 
process generated by A concentrates on the absorbing set stf and has ini- 
tial condition given by the harmonic measure v x * if X K (0) — > x* /or some 
x* E E. 

Proof. Convergence of the semigroups to the semigroup with generator A 
follows from [14J, Theorem 2.1. The only condition to check is dom(A) C 
&(\ — A) for some A > 0. This condition is clearly satisfied if A generates 
a Markov process, because in that case M{\ — A) equals C(&f,W). 
With Dyson's formula ([T7j p. 367) we have for all / in the core 

S K (t) f - e KtL ' f = f e^ L 'LS K {s)fds . (47) 
J o 

Using convergence of the semigroups S K (t) and e KtL ' as k — > oo we get 

S(t)f-Pf= f AS(s)fds . (48) 
J o 

This shows that A generates S(t) and in particular that S(0) = P. For har- 
monic / E J^f(E) we have S(0)f = Pf = f and S(t) is indeed a semigroup 
on J(f(E). This implies that the limit process generated by A concentrates 
on the absorbing set si as explained in more detail below, and that by the 
definition of the projection P (|45[) the initial condition of the process is given 
by the harmonic measure v x . □ 

Every L '-harmonic function on E is uniquely determined by its boundary 
values on srf (cf. (]4"5])). So in addition to the projection P we can define the 
harmonic extension 

& : C(^,R) -> 31? (E) where &f(x) := / v x (da)f(a) . (49) 
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Notice that P and are denned by the same formula, but act on different 
spaces. In particular we have Pf = (=^(/|^)), i-e., the harmonic projection 
is the same as the harmonic extension of the restriction of / to s/. With 
this notation we can define 

A' : C(sf, R) -»• C(s/, M) where A' f := {A3? f) \* , (50) 

which highlights the fact that the limit semigroup generated by A actually 
defines a process on the absorbing set s/. Via the identification of J^{E) 
with C(s/, R) the operators A and A' both describe infinitesimal motion 
away from s/ due to L, instantaneously followed by projection onto s& ' . 
Formally, in order to define A 1 on / € C(A, R) and act with L we first 
have to extend the function / to the whole configuration space E: therefore 
we first harmonically extend / (action of then apply L, then project 
(action of P) and restrict again to s/ ' . 

In order to apply Theorem 14.11 we have to assure that the operator A' 
generates a Markov process. This will be done in the next Section and below, 
where we explicitly compute A 1 and recognize it as the generator of a process 
on s/ with a jump and a diffusion part. The condition of convergence of the 
semigroup of the fast diffusive part of the generator in our case is covered 
by the following result. 

Proposition 4.1. Let L' = | Y^i^sPih^yiVj^j be the generator ^ of 
the Wright-Fisher diffusion X' with absorbing states s/ ([5^. Then 

\\e tL ' f - P f\\oo ^0 ast^oo. (51) 
Proof. Using that P is the projection on the absorbing set we have 

{e tL> - P)f(x) = E x [f(X'(t)) - /(X'(oo))] • 
Denoting by r the time of absorption of the Wright-Fisher diffusion, we get 

E x [f(X'(t)) - /(X'(oo))] = E x [(f(X'(t))-f(X'(r)))l(r>t)] 

< 2||/|| 00 P ie [T>t]. (52) 

This converges to zero since / is bounded on the compact state space E, 
and the Wright-Fisher diffusion gets absorbed almost surely in finite time 
with a uniformly bounded mean absorption time (see e.g. [21]). □ 

To illustrate how the projected generator A' of (|50p looks like in a con- 
crete case, let us look back at the example of two sites with E = [0, 1] , 
studied in the previous section. The X'-process is the Wright-Fisher diffu- 
sion with generator L 1 = 2x(l — x)d x , with absorbing states and 1. The 
X process is the deterministic process with generator L = —^{2x — l)d x . 
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Since X'{t) is a martingale, stopping at the time of absorption gives that 
the harmonic measure equals 

v x {{\}) = 1 - v x ({0}) = x . 

As a consequence, for a function / : {0, 1} — > R, its harmonic extension 
is simply given by linear interpolation: 0>f(x) = xf(l) + (1 — x)f(0) and 
therefore 

(L^f)(x) = -^^^(*/(l) + (l-x)/(0)) = -^^(/(^-/(O)) • 
Thus 

Af(x) = (PL0>)f(x) = - [x(f(0) - /(l)) + (1 - z)(/(l) - /(0))] 
which takes values 

= /(0) - /(l) and A/(0) = /(l) - /(0) . 

So A is the generator of the process which flips at rate 1 between the states of 
the absorbing set srf = {0, 1}, as we found in Section T3.2I using a martingale 
characterization. 



5 Computation of the generator A: the case of 
three sites 

In the case of general connections p(i, j) and general finite sets S, the possible 
limiting behaviors of the process will become more complicated. Therefore 
we focus on the case of Theorem 12.21 with p(i,j) € {0, 1} and compute the 
generator A (f50|) using that 

fe hL -I\ fe hL -I\ 

PL = P lim = lim P (53) 

h\o \ h ) h\o \ h ) 

by continuity of the projection P. We start here with a system of three sites 
with closed boundaries and nearest neighbor connections. This case already 
contains most of the interesting aspects of the limiting dynamics which now 
has both a diffusive and a pure jump part. 

5.1 Generator and absorbing set 

We consider the three site case with p(l, 2) = p(2, 1) = p(2, 3) = p(3, 2) = 1 
and p(l, 3) = p(S, 1) = 0. The generator ([6]) is 

Ln = -a(x 1 - x 2 )d 12 - a(x 2 - x 3 )d 2 3 

+ -^x x x 2 d\ 2 + — x 2 x 3 d 23 (54) 
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where we write a = a/4 in the following to simplify notation. Recall the 
notation from © for this case: the "fast Wright-Fisher diffusion" has gen- 
erator 

L' = \ Xl x 2 d\ 2 + \x 2 x,dl, , (55) 

whereas, the "slow deterministic" part has generator 

L = -a(xi - x 2 )di2 - a(x 2 - x 3 )d 2 3 ■ (56) 

The process x'(t) with generator V is a Wright-Fisher diffusion, for which, 
x±(t), x' 2 (t), x' 3 (t) and Xi(t)x' s (t) are martingales since L'x\ = L'x 2 = L'x^ = 
L'{x\Xz) = 0. Its absorbing set is 

n/ = (0,z,0) U {(u,0,w) : u + w = z) (57) 

given by a point and a line segment. Here we assume that the total mass 
in the system is z < 1, for our results to generalize more easily to larger 
systems. 

5.2 Dynamics on the line segment: diffusion and jumps 

We start with a point on the line segment x(0) = (u,0,w) and need to 
compute (see ([35]) ): 

Af(x(0)) = lim 1 (E xW [f(x'(r))} - f(u, 0, w)) (58) 

where x(h) = (u(l — ah), zah, w(l — a/i)) € E \ stf is the initial condition for 
the Wright-Fisher diffusion x'(t). This initial condition comes from starting 
from the point (u, 0, w) from the absorbing set and evolving according to 
the deterministic part of the dynamics a small amount of time h. 
Then, since x' 2 (t) € {0, z} takes only two values, we have 

E x{h) [f(x'(r))) = f(0,z,0)F x{h) [x' 2 (T) = z] 

+E x{h) [f(x'(T))l(x' 2 (T) = 0)} . (59) 

Since x' 2 (t) is a martingale on [0, z], the probability in the first line is given 
by ah and we have 

E*(h) [f(x'(r))] - f(u, 0, w) = ah(f(0, z, 0) - /(«, 0, w)) 



+E 



x(h) 



{f(x'(r)) - f(u, 0, w)) l(4(r) = 0) . (60) 



The first line corresponds to a jump part, where all the mass z ends up in 
site 2. The second line has to be a continuous part, since the probability of 
the event x' 2 (t) = is of order 1 + o(h). Since Xi(t) is a martingale we have 

E x(h) Wi(t)-u] = -uah = -uah + [(x[{t)-u) 1{x' 2 {t) = 0)] , 
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where the second equality uses again ¥ x ^[x' 2 (t) = z] = ah and the fact 
that x' 2 (t) = z implies x'i(r) = 0. Therefore 

E x(h) [{x\ (r) - u) l(x' 2 (r) = 0)] = (61) 

and similarly for x' 3 (t), and the drift of the continuous part vanishes. It 
remains to compute the diffusion part. 

E s(h) [(xi(r) - u) 2 ] = E x(h) [x'^rf] - 2uE x(h) [x'^r)] + u 2 

= E x{h) [x[(t)(z - x' 2 (r) - 4(r))] - 2n 2 (l - ah) + u 2 

= uz(l - ah) - - uw(l - a/i) 2 - 2u 2 (l - ah) + u 2 

= ahuz + 0{h 2 ) , (62) 

where we have used that x / 1 (t)x2(t) = and that x'i(t)x' 3 (t) is a martingale. 
Again splitting the expectation with respect to the value of x' 2 (t) we get 
analogously to the above 

E x i h \ [(x' x (r) — u) 2 1(x 2 (t) = 0)] = ahuz — ahu 2 = ahuw . (63) 

The same holds for x' 3 (t), and obviously E x q^ [x 2 (T) k 1(x' 2 (t) =0)] =0 for 
all k = 1,2 To get the covariances we compute 

Ex(h) [(x' 1 (r)-n)(x / 3 (r)- U ;) 1(4 (r) = 0)] = 

= E^) [xi(t)4(t)) 1(4(t) = 0)] - ««; = -2a to , (64) 

where we used (|61|) and that x'^t) x' 3 (t) is a martingale. Covariances with 
x' 2 (t) again vanish. 

Now, in (|6TI|) we have up to second order 

3 

f(x'(r))-f(u,0,w) = Y / d x J(x(0)){x l (r)-x l (0)) 

i=l 

1 3 

+ 2 E ^ l ^,/(x(0))(x i (r) - x f (0)) (x,(r) - ^(0)) , (65) 

where higher order terms are of order o(h) and not relevant. Taking expec- 
tations, the first order terms vanish, and the second order diffusive terms 
can be written as (recall that a = a/4) 

E x(h) [{f(x'(r)) - f(u, 0, w)) 1(x' 2 (t) = 0)] = ^ahuwd 2 3 f(u, 0, v) + o(h) 

using the same notation as in ([5]). Plugging this and ()60p in (|58p we obtain 
for the limiting generator A: 

Af(u, 0, w) = a(/(0, z, 0) - f(u, 0, w)) + ^auwd 2 n f(u, 0, w) , (66) 

consisting of a jump part and a Wright-Fisher diffusive part with effective 
diffusivity auw. Notice that as soon as u or w = the diffusive part vanishes, 
and the generator consists purely of the jump part where the total mass z 
moves from site 1 or 3 onto site 2. 
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5.3 Dynamics of a single condensate: jump process 

For the other possible initial condition x(0) = (0,z,0) we have x(h) = 
(zah, z(l — 2ah),zah)} as initial condition for the Wright-Fisher diffusion. 
Again we have 

E x{h) [f(x'(r))] = f(0,z,0)F x{h) [x' 2 (r) = z] 

+E x{h) [f(x'(r))i(x' 2 (r) = 0)] , (67) 

where this time the probability in the first line is of order 1 — 2ah. Since 
x' t (t) X3 (t) is a martingale we have 

E x{h) [x[(T)x' 3 (T)] ={zahf (68) 

which implies that the probability to end up in a state (u, 0, w) with u, w > 
is negligible. Therefore 

W{t)) l(4(r) = 0)] = 2ah(±f(z, 0, 0) + ^/(0, 0, z)) , (69) 

and the generator consists only of jump parts which can be written as 

Af(0, z, 0) = a{f(z, 0, 0) - /(0, z, 0)) + o(/(0, 0, z) - /(0, z, 0)) . (70) 

Therefore, once all the mass concentrates on a single site this remains, and 
the corner points ^ are absorbing for the limit process. 

The closure of the operator A defined by (fT0"|) . (f66"|) clearly is the generator 
of a well-defined Markov process on the set g/ defined in (|57p . Furthermore, 
tightness of the sequence of processes generated by Ln ([54"]) is proved in 
Lemma 14.21 Therefore, we can conclude for the three site case, that the 
processes generated by (f54"|) indeed converge to the process generated by the 
projected generator A defined by (fT0|h (f66|h This concludes the proof of 
Theorems 12.11 and 12.21 for the case of three sites. 

6 The general case continued 
6.1 Motion of a single condensate 

First we show that once a point € ^ is hit in the limit process, i.e., the 
situation physically corresponding to a single condensate, then the process 
becomes a pure jump process on the corner set ^ = {e^ : i £ S}. This in 
combination with Lemma 14.11 finishes the proof of the first item in Theorem 
12.11 Recall the decomposition = L + 9j^L' of the generator of the 
auxiliary process given in ([6]). 

Lemma 6.1. 1. For every Borel measurable subset K C E\ c € we have 
for the function v.{K) ■ x 1— >• v x {K) 

(LvXK))( ei ) = 0. (71) 



24 



As a consequence, Af(x) = for all x G E\ c € and hence, 



Af(e i ) = Y,^P(iJ)(f(e J )-f(e l )) 



i.e., the projected process is a pure jump process on ^ . 

Proof. Fix K C E\% '. If p(i,j) > then starting from a point x & E 
such that Xk = 0, k {i,j}, the Wright-Fisher diffusion L' remains in the 
plane Eij := {rei + sej : r, s £ R} and hence is absorbed at or ej. As a 
consequence 



Since K was an arbitrary Borel measurable subset of E \ we conclude 
for every bounded Borel measurable function / with support contained in 
E\^ that 



V Js2 / 
Further, if i {£;, /} then obviously, for any smooth function / : — > M 



because both and are zero when evaluated at x = ej. 

Therefore, since L in Q is a linear combination of terms of the form 
~j( x k — xi)dki, we conclude that for any smooth function / : — >■ M with 
support contained in E \ ^€ 



which implies 



(dijv{K))(ei) = . 






Therefore, we obtain for all / G C(si ,IR) 




(72) 



since f v v ei (da) f (a) = /(ej). □ 
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6.2 The general case of Theorem 12.21 



The proof of the general case is a combination of the arguments used in the 
three site case discussed above, and we will be brief in this section avoiding 
a full treatment with a lot of technical notation. Recall that p(i,j) £ {0, 1} 
and let x(0) € £/ be a general initial condition, for which Xi (0)^(0) = 
whenever p(i,j) = p(j,i) = 1- Again, under the Wright-Fisher process 
(x'(t) : t > 0), with generator V = (1/2) EijP()i)K%^> a11 X 'M and a11 
products a^(i)a^(i) for which p(i,j) = are martingales. Analogously to 
the computations in Section after the action of the deterministic drift 
generator L = ~ a P(h j)(Ui ~ Vj)^ij f° r infinitesimal time h we have 
an initial condition x{h) for the diffusion L', and denote by r the time of 
absorption. 

The first observation is that at most one initially empty site can gain a 
macroscopic amount of mass through a jump. 

Lemma 6.2. Xj(0) = and x^(r) > is possible for at most one i € S, and 
in that case YljeS p(*> j) x 'j( T ) = 0- 

Proof. Let i and j be two initially empty sites for which Xi(h), Xj(h) = 
O(h). Then, if p(i,j) = 1 we have x^(r)x^(r) = since r is the hitting time 
of the absorbing set so that at most one of them can take a positive 
value. Similarly, if p(i,j) = we have x^(r)x^(r) = 0(h 2 ) since the product 
is a martingale, and since we take the limit h — > again at most one can 
be positive after the diffusion is absorbed. Since the pair i,j was arbitrary, 
this implies that under the action of the generator A at most one empty site 
in the system can gain mass. Then at the time of absorption in all other 
connected sites have to be empty. □ 

Secondly, jumps of a macroscopic amount of mass occur only onto neigh- 
bouring empty sites. 

Lemma 6.3. If x0) = and J2jesP(^j) x j(^) = °> ^ en x 'i( T ) = °- V 
Xi(0) > then x[{t) < Xi (0) + 0(h). 

Proof. By the condition site i has no direct neighbour with xj(0) > 0, 
and therefore if Xj(0) = under the action of generator L in time h we get 
Xi(h) = 0(h 2 ). Then since x[{t) is a martingale, its mass remains negligible 
after the Wright-Fisher diffusion V gets absorbed and x-(r) = 2j(0) = 0. 
If Xj(0) > we also have X^esP^j j) x j(®) = and thus Xi(h) < Xj(0) + 
0(h 2 ). Then the martingale property of x[ implies that absorption at 
x'^t) = Xi(0) + 0(1) has negligible probability of order h 2 . □ 

Thirdly, mass between two macroscopically occupied sites moves at most 
continuously. 
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Lemma 6.4. If Xi(0)xj(0) > 0, then either sj(r)x'-(r) = or x-(r)x^(r) = 
a* (0)^(0) + 0(/i). 

Proof. Xi(Q)xj(0) > implies p(i,j) = 0. In the first case the mass of one 
of the sites jumped to another site. Now assume this is not the case and 
note that x-(i), x'j(t) and x^ifyx'jit) are martingales. Let f,g : [0, 1] — > R be 
two smooth functions, then 

L '{f(xi)g(xj)) = ^2 (p(i,k)xix k df k f(xi) + p(j, k)xjX k d% k g(xj)J . (73) 
kes 

Then Xi(h)xj(h) = O(l) which implies that x k (h) = 0{h) for all k connected 
to either i or j and 

L'(f(x i (h))g(x j (h)))=0(h) . (74) 

So either a neighbour x' k (t) gains a macroscopic amount of mass during the 
Wright-Fisher diffusion, which results in the first case of one of the piles 
jumping to k, or we have 

EswiMWWjir))] = fixiiOMxjQ)) + 0(h) . (75) 

Since / and g are arbitrary, the latter case implies the second statement. 
□ 

Lemma 6.5. If Xj(0) = and x[(t) > then x-(r) = J2jesP(j^) x j( ) + 
O(h). 

Proof. This follows from Lemmas 16.21 and 16.41 i is the only site in the 
system where a macroscopic amount of mass jumped to and its neighbours 
are empty at absorption, therefore it has to absorb all the original mass of 
its neighbours up to amounts of order h which have been shared with other 
sites. □ 

This characterizes all the possible jump events under the generator A. 
We see that mass can only jump onto a single empty site j which collects all 
the surrounding mass, and the rate for this event can be computed analo- 
gously to the previous section to be j Y2iesP(h j) x i(ty- ^ n parallel to jumps, 
mass can move continuously between macroscopically occupied sites, which 
we describe next. We will use the notation p(i,j) introduced in (fT5|) for 
the number of two-step connections between two sites i,j 6 S. The crucial 
ingredient to conclude is additivity of the jump rates of the Wright-Fisher 
diffusion in each component, as is described below. 

Lemma 6.6. Assume that p(i, j)xi(0)xj(0) > then the generator A con- 
tains a term ap{i,j)xi{fS)xj(^S)d^A and the continuous part of A is given by 
the sum of all such contributions. 
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Proof. For i,j 6 S as given, assume that x^(r)^(r) > so that none of 
the two piles jumps and the total mass is conserved up to order h, which 
has 0(1) probability. Then the system on sites i, k,j with intermediate site 
k can be seen as an effective three site system as studied in the previous 
section. If there is a unique intermediate site, the diffusive part of the 
generator follows analogously to this case. If there are several intermediate 
sites k, the different paths between i and j can be viewed as independent, 
since the diffusion rates of the Wright-Fisher generator acting on sites i and 
k given by 

L'ik = XiX k df k (76) 

are proportional to x k (similarly between j and k). Therefore the total 
effective Wright-Fisher diffusion rate between sites i and j is simply given 
by a sum over all connections, as given by p(i,j) in f j 1 3 p . Since the above 
rates are also linear in Xi (and Xj, resp.), the same additivity applies if i 
or j are connected to more different sites via two-step connections. This 
leads to independent contributions for each such connection as claimed in 
the Lemma, which are summed over to give the full continuous part of the 
generator. 

It remains to show that there is no mass exchange over distances more 
than 2 steps. This is most easily seen in a system with 4 sites S = {1, 2, 3, 4} 
with nearest neighbour jumps and initial condition (u, 0, 0, w) and u+w = z. 
Analogous to the computation in ([62]) for three sites we get 

^x(h) [(^'i( T ) - u ) 2 ] = ah(uz - uw) + 0{h 2 ) . 
This implies analogously to (JMJ) that 

E i(/i) [(i' 1 (t) - n) 2 l(4(r) = 0)] = ah(uz - uw) - ahu 2 + 0(h 2 ) = 0(h 2 ) , 

so the change of x\ is not on a diffusive scale and negligible. The same holds 
for X4, and this argument can be easily generalized to arbitrary S. □ 
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